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Massive multiple-input multiple-output (MIMO) is becoming a key 
technology for future 5G wireless communications. Channel feedback 
for massive MIMO is challenging due to the substantially increased 
dimension of MIMO channel matrix. In this letter, we propose a 
compressive sensing (CS) based differential channel feedback scheme 
to reduce the feedback overhead. Specifically, the temporal correlation 
of time-varying channels is exploited to generate the differential channel 
impulse response (CIR) between two CIRs in neighboring time slots, 
which enjoys a much stronger sparsity than the original sparse CIRs. 
Thus, the base station can recover the differential CIR from the 
highly compressed differential CIR under the framework of CS theory. 
Simulations show that the proposed scheme reduces the feedback 
overhead by about 20% compared with the direct CS-based scheme. 


Introduction: Massive multiple-input multiple-output (MIMO) is 
emerging as a key technology for future 5G wireless communications, 
ly-^por massive MIMO systems, accurate channel state information at 
T __j:ransmitter (CSIT) is essential to improve both the spectrum and 
f~y nergy efficiency (T). One approach to acquire CSIT is to adopt time 
^odivision duplexing (TDD) mode by exploiting the channel reciprocity. 
However, as frequency division duplexing (FDD) outperforms TDD in 
—jielay-sensitive and traffic-symmetric applications, FDD still dominates 
I cellular networks without channel reciprocity f2j, where effective 

channel feedback is crucial to enable the expected system performance of 
^“'massive MIMO. As the dimension of the MIMO channel matrix becomes 
Very large in massive MIMO systems, while the dedicated resource for 
ch annel feedback is limited in the uplink, channel feedback becomes a 
challenging problem. Conventional channel feedback methods based on 
p odehook are not suitable for massive MIMO systems, since the size of 
the codebook has to be extensively expanded to guarantee an acceptable 
J3SIT accuracy (3). For non-codebook based channel feedback, which 


Utsuallly requires large overhead for channel feedback, compressive 
sensing (CS) has been applied to reduce the feedback overhead in 
massive MIMO systems o-dj. Some CS-based channel feedback 
^schemes (3], (4) exploited the spatial correlations among antennas to 
OOtompress the MIMO channel matrix for overhead reduction, but MIMO 
^—Channels are not always spatially correlated in practice, especially when 
\Qi istributed MIMO is considered GQ. In a, the channel impulse response 
T t^j"(CIR) is directly compressed by exploiting the time-domain channel 
C ~ 'S parsity to reduce the feedback overhead, but the overhead is still high, 


r- 


Especially when the CIR is not sufficiently sparse. 


In this letter, we propose a CS-based differential feedback scheme 


i"TTto reduce the feedback overhead by exploiting the temporal correlation 
of MIMO channels. Specifically, as the temporal channel correlation 
. .exists in both centralized and distributed MIMO systems, we propose to 
^generate the differential CIR by differentiating the CIRs in two adjacent 
^ime slots, which enjoys a much stronger sparsity than the original 
CIRs. We then propose to compress the differential CIR with lower 
$—compression ratio, which equivalently requires much lower feedback 
^Mtverhead. Simulation results show that to achieve the same normalized 
mean square error (NMSE) performance, the proposed scheme can 
reduce the feedback overhead by about 20% than the direct CS-based 
channel feedback. 


Temporal correlation of time-varying MIMO channels: In the t-th time 
slot, the CIR between the n-th transmit antenna at the base station (BS) 
and the single-antenna user in massive MIMO systems can be denoted 
as = \hn\o), h„\ 1), • ■ ■ , h$(L — 1)] for 1 < n < N, where N 
is the number of transmit antenna and L is the maximal channel delay 
spread. The CIR is often considered to be sparse (i.e., there are K 
non-zero elements in , 1 < K -C L) since it usually contains a few 
significant paths in broadband communications (6j. The sparse CIR series 
{h^ }J = i in T consecutive time slots exhibits temporal correlation even 
when the MIMO channels are fast time-varying 0. We model such 
temporal correlation through a variation of the CIR’s support (i.e,, the 
positions of non-zero elements) and an evolution of the amplitudes of 


non-zero elements. To be specific, the time-varying sparse CIR can be 
characterized by a support vector pit' 1 and an amplitude vector as 
follows 


h (t) - n (t) o a (t) 
ll n — Pn u “n j 


(i) 


where p^ (I) € {0,1}, a ^ ( l ) G C and o denotes the Hadamard product. 
To model the variations of the CIR’s support over time, the 

Z-th elements of T support vectors {Pn^(()}t=i ' n consecutive T 
time slots can be modeled as a first-order Markov process m 
which can be fully characterized by two transition probabilities pro = 
Pr{Pn +1) (0 = l|Pn ) (0 = 0} and poi =Pr{p£ +1) (0 = Ob^Z) = 1}, 
and a distribution = Pr{p^(Z) = 1} in the initial time slot t=l. 
For the steady-state Markov process, where Pr{p^(Z) = 1} = p,Vt,n, 
only two parameters poi and p are enough to characterize such process 
since pro = ppoi/(l — p). Then, to model the evolution of the CIR’s 
amplitudes over time, the amplitude vector a^ can be modeled by the 
first-order autoregressive model as follows 0 

= pal t_1) + \/l - p 2 w (t) , (2) 

where p= is the correlation coefficient given by the zero- 

order Bessel function of the first kind with /^ being the maximal Doppler 
frequency and r being the time slot duration, the parameter w W denotes 
a noise vector independent of a£ 4 1 \ and its entities are independent 
and identically distributed (i.i.d.) zero-mean complex Gaussian variables 
following the distribution CAf(0, o^). 


Compressive sensing based differential channel feedback: As shown 
in Fig. 1, unlike the direct CS-based channel feedback scheme, which 
directly compresses the sparse CIR with a sensing matrix, the proposed 
CS-based differential channel feedback scheme further exploits the 
temporal correlation of MIMO channels to generate the differential CIR 
which enjoys a much stronger sparsity than the original CIRs. The 
differential CIR is then compressed by the sensing matrix, which is 
known to both the BS and the users, and then CS algorithms can be 
used to produce a precise recovery of the differential CIR after channel 
feedback. In order to avoid feedback error propagation, we further 
propose that users execute initialization after a certain (either fixed or 
adaptive) number of time slots, whereby the original CIR is compressed 
with a relatively higher compression ratio to ensure a precise initial CIR 
recovery. More specifically, the procedure of the proposed scheme can be 
decomposed into the following three steps. 


User 


BS 



Fig. 1. Compressive sensing based differential channel feedback. 

1. Differential operation: Unlike conventional feedback schemes which 
directly feedback the estimated channel in the current time slot, the 
proposed scheme aims to create a more sparse representation of the 
channel by computing the difference between two estimated CIRs 
hn 11 and in the previous (f — l)-th and the current f-th time 
slots, respectively. Mathematically, the differential CIR Ah^’ can be 
calculated as 

Ah^hW-hf 1 ’ 

= pi f) ° (a? - _1) ) + (p n' 1 - p£ t_1) ) ° (3) 

= Pn ) ° [ 1 /l-p 2 w-(l-p)a^' 1) ] + (p^ ) -pL t_1) ) °an _1) - 

The first item on the last line of (3) is vety small due to the smooth 
evolution of the CIR’s amplitudes, e.g., pta 1 (8). The number of non¬ 
zero elements of the second item is also small (e.g., not larger than K ) 
since the sudden disappearance and appearance of a significant path do 
not occur often. Fig. 2 presents the snapshot of the previous CIR, the 
current CIR, and the differential CIR of the channels described by (1) 
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and (3) with the following parameters: L = 200, poi =0.05, /r = 0.1, 
fd = 10Hz, t = 1ms, and a w = 1. The initial amplitude is a random 
vector, the entries of which are i.i.d. complex Gaussian variables with 
zero mean and unit variance, and the initial support p,^ is a random 
vector, the entries of which are i.i.d. Bernoulli variables. We can observe 
from Fig. 2 that the differential CIR enjoys a much stronger sparsity 
than the original CIR due to the temporal correlation of channels, which 
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Fig. 2 A snapshot of the previous CIR, the current CIR, and the differential 
CIR: (a) the previous CIR: (b) the current CIR: (c) the differential CIR. 

2. CIR compression and recovery: CS theory shows that a sparse signal 

Ah^ can be compressed by a sensing matrix # w j t [ 1 ^ ^ 

L into a low-dimensional measurement vector y = •f’AhJ,*' (y will be 
fed back to the BS). If <1? satisfies the restricted isometry property (RIP) 
|91 , the sparse signal Ah^ can be reliably recovered by the BS from 
the noisy measurements y = Ahj,^ + n. where n denotes the noise 
of feedback channel, the entities of which are i.i.d. complex Gaussian 
variables following the distribution CJ\f(0, a^f). In this letter, we set •!? as 
the commonly used random Gaussian matrix da. As for the CS signal 
recovery algorithm, we adopt the widely used algorithm called subspace 
pursuit (SP) due to its low complexity and robustness to noise oa. 

3. Error propagation avoidance: Precise recovery of the CIR in the 
initial time slot is of great importance since the channel feedback 
error in early stage will propagate in the following CIR recovery 
process. Additionally, an unexpected recovery error of the differential 
CIR will also impair the subsequent CIR recovery process. To avoid 
such error propagation, the proposed feedback scheme will be reinitiated 
every P time slots. Although such initialization will cause relatively 
large feedback overhead to ensure a precise recovery of the initial 
CIR, it will be verified numerically in the next part that the highly 
compressed differential CIR will save much more feedback overhead 
in the subsequent time slots. Thus, the total feedback overhead of the 
proposed scheme is still lower than the direct CS-based feedback scheme. 

Simulation results: Simulations have been conducted to validate the 
NMSE performance of the proposed CS-based differential channel 
feedback scheme. We consider a massive MIMO system with the same 
parameters used for the generation of Fig. 2 with N = 32 antennas. The 
initial amplitude vector aj, 1 ' and support vector p,^ for 1 < n < N are 
independent of each other. We denote the channel compression ratio (i.e., 
relative feedback overhead) as p = M/L, and two cases are considered: 
p = 25% and p = 45%. For the proposed scheme, p is denoted as the 
average compression ratio over P feedback slots. The compression ratios 
in the initial time slot and subsequent time slots are 45% and 15% for 
the first case, and 65% and 35% for the second case. Thus, the average 
compression ratio of the proposed scheme is 25% for the first case 
and 45% for the second case when P = 3 is considered as an example 
(Note that P can be adaptive to the channel condition). As shown in 
Fig. 3, the proposed scheme outperforms the direct CS-based scheme 
when the same compression ratio is considered in both cases. In the 
first case, the direct CS-based scheme fails to work since the feedback 


overhead is insufficient, while the proposed scheme performs well. In the 
second case, the proposed scheme achieves a 4 dB SNR gain compared 
with the direct CS-based scheme. Additionally, the proposed scheme 
with p = 25% and the direct CS-based scheme with p = 45% have very 
similar NMSE performance, which indicates that the proposed scheme 
can significantly reduce the feedback overhead from 45% to 25% to 
ensure a reliable channel feedback. 



Fig. 3 NMSE performance comparison between the direct CS-based feedback 
scheme and the proposed feedback scheme. 

Conclusion: This letter investigates the challenging problem of channel 
feedback in massive MIMO systems. It is found that by exploiting the 
temporal correlation of time-varying channels, the generated differential 
CIR between CIRs in neighboring time slots enjoys a much stronger 
sparsity than original CIRs. We have shown that the proposed CS- 
based differential channel feedback scheme can achieve better NMSE 
performance than the direct CS-based channel feedback scheme, and 
reduce the feedback overhead by about 20% for massive MIMO systems. 
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